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In this chapter we review the present status of the Gaussianity studies of the CMB anisotropics, in- 
cluding physical effects producing non- Gaussianity, methods to test it and observational constraints. 
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I. INTRODUCTION 



The study of the CMB temperature and polarization anisotropy has had an essential role in establishing the standard 
cosmological model: a flat ACDM model composed of baryonic matter (« 4%), dark matter (« 24%) and dark energy 
(« 72%), with nearly scale-invariant adiabatic fluctuations in the energy-density. The information extracted from the 
CMB to derive the cosmological model has been mostly based on the CMB anisotropy power spectrum. However, much 
more information is still waiting to be extracted from the statistics of the CMB beyond the second-order moment. 



* Chapter of the book "Cosmic Microwave Background: from quantum fluctuations to the present universe" , IAC Winter School 2007. 
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The study of the Gaussianity of the CMB anisotropics has become in recent years a very relevant topic in the 
CMB field. The main reason for that is the availability of high sensitivity and high resolution maps provided by 
the new generations of CMB experiments, with which the Gaussian character of the anisotropics predicted by the 
standard inflationary paradigm can be tested at a high accuracy. CMB anisotropics form an ideal data set for 
testing the Gaussianity of the primordial energy-density fluctuations since the dominant physical effects producing 
the anisotropics involves just linear physics, and Gaussianity is preserved under linear transformations. Another 
approach to test it is by studying the large scale structure of the universe as described by the galaxy distribution. The 
density perturbations in the linear regime should also be a good representation of the initial conditions. However, 
galaxy formation is a very non-linear process involving complex physical effects which very much complicate the 
analysis. 

Different models of the early universe have been proposed to naturally account for the early stages of its history 
and, thus, not to have to rely on ad hoc initial conditions. They are based on different theories, like string theory or M 
theory, and include or not an inflationary phase (an example of the latter being the Ekpyrotic cosmology [HI]). Most 
of those models predict very specific properties about the probability distribution of the CMB anisotropics, which in 
many cases imply deviations from Gaussianity with amplitudes within reach of present or near future experiments. 
Some of those models are already constrained by present data and many more are expected to be disproved in the 
coming years, specially with the launch of the Planck satellite at the end of 2008 (or beginning of 2009). 

The analysis of the angular distribution of the CMB anisotropies is a formidable task with profound consequences 
on our understanding of the universe. Before any meaninful result can be achieved, it is, however, crucial to control 
at a high level of precision all possible systematics that can be introduced by the experiment and the pipeline process 
used to reduce the data. On the other hand, there are several Galactic and extragalactic emissions in the microwave 
band which blur the CMB signal. Disentangling the cosmic signal from the others is also crucial and requires a very 
good knowledge of the astrophysical emissions. A lot of observational and theoretical effort has been dedicated to 
that aim in recent years. Below we summarize the present status of the most important aspects of the analysis of 
Gaussianity and the results achieved. 



II. THE ISOTROPIC GAUSSIAN RANDOM FIELD (IGRF) 



One of the most robust predictions of the standard inflationary model is that the CMB anisotropies should be well 
represented by an isotropic Gaussian random field (IGRF) on the celestial sphere. This is a very powerful prediction 
to test the standard model. There are many statistical quantities for a IGRF that can be derived analytically, which 
greatly facilitates the Gaussianity test of the CMB data. However, a disadvantage may come from the central limit 
theorem that implies that the sum of several independent non-Gaussian distributions tends to a more Gaussian one. 
As we will see below, this complicates the Gaussianity analysis since the observational data is composed of several 
contributions which can be either intrinsic or extrinsic to the CMB anisotropies. 



A. Definition 



A random field defined on a given support is said to be Gaussian if for any N points of the support x\, ...,x„ the 
values of the random field yi, ...,y n follow a multinomial distribution 

/(fi'-'Vn) = ( 27r) »/2| M |i/2 ex P < W >) M ^iVo~ < Vi >)j C 1 ) 

where M is the covariance matrix defined as My =< (j/j— < j/j >)(l/j~ < Uj >) >■ Thus the n-pdf (probability 
density function), and also all the moments, is given in terms of just the first two moments. 

In the case of the temperature anisotropies of the CMB the mean value is set to zero, and thus the standard model 
predicts a Gaussian pdf characterized by only the second moment. The support is in this case the 2D sphere. Thus 
the temperature anisotropies can be expanded in terms of the spherical harmonic coefficients, 

AT 

A(n) = — (nJ^aAW (2) 

Ira 

where Y/ m (n) are the spherical harmonic functions for direction n and ai m are the spherical harmonic coefficients 
which for the standard model are Gaussian distributed . 
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If the field is isotropic then the two-point correlation function only depends on the modulus of the difference of the 
two directions. In harmonics space, isotropy translates in that the harmonic coefficients are uncorrelated 

< O'lmO'l'm' >= CiSu'S mm ' (3) 

where Sij is the Kronecker delta. Gi is the temperature power spectrum which for a realization can be estimated as 

B. Properties 

The IGRF is one of the best studied random fields and many of its properties have been thoroughly analysed. One 
of the most remarkable properties is that the expectation of any even combination of the field A(n) can be given in 
terms of the second moment, the 2-point correlation function, and the expectation of any odd combination is zero; i.e. 
the n-point correlation functions for n odd are zero. The same property translates to the spherical harmonics space 
where the expectation of even combinations of the coefficients ai m can be expressed in terms of the power spectrum 
C[ and the expectation of odd combinations, e.g. the bispectrum, is zero. More generally, a very useful characteristic 
of a IGRF is that the expectations of many statistical quantities can be calculated (semi) analytically. This is the case 
for the morphological descriptors, number, shape and correlation of peaks (maxima), scalars on the sphere, ... 

The Minkowski functionals are useful descript ors of the morphology of point sets or smoothed fields in spaces with 
arbitrary dimension d (see e.g. [l|). As stated in [138], under a few simple requirements any morphological descriptor 
can be written as a linear combination of d + 1 Minkowski functionals. For the sphere they are therefore 3; namely, 
total length of the contour C(v) of the excursion set, total area A{v) of the excursion set and the genus G(v) above 
a threshold v. Their average value per unit area for a IGRF can be simply given by 

<W>=iexp(4) , (4) 



1 2 



< A(v) >= ~ - -j= / exp (-x 2 ) dx , (5) 



(6) 



where 6 C = ( — ^tt^jj) 1 ^ 2 is the coherence angle of the random field which is defined by the ratio of the two-point 
correlation function C(8) to its second derivative at zero lag (see Fig. [1]). 

Properties of peaks in a 2D IGRF have been studied by [llOj ] , [23| and Q . The local description of a peak involves the 
second derivative of the field along the two principal directions. The curvature radii are defined in the usual way from 
the second derivative of the temperature anisotropies A at the position of the maximum: R\ = [— A" (max)/ 2a q] ~ 1 / 2 
i?2 = [— A2 , (maa;)/2(To]~ 1 / 2 , where ctq is the anisotropy rms and A" is the second derivative along the principal 
direction i. The two invariant quantities, Gaussian curvature k and eccentricity e can be constructed from them: 



R\R-2 



Ri 



1/2 

(7) 



It is straightforward to obtain the number of peaks on the celestial sphere N(k, e, v) dn de dv with Gaussian curvature, 
eccentricity and threshold between (re, k + dn), (e, e + de) and {y, v + dv), in terms of the two spectral parameters 7 
and that characterize the cosmological model: 

7 = ^ , 0* = 2V^i , (8) 
a 2 = C(0) , <rl = -2C"(0) , «7 2 2 = |(7W(0) . (9) 
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FIG. 1: Expected value of the Minkowski functionals for an IGRF obtained from eq.|4] eq. eq.[6] Also plotted are the average 
values (asterix) and error bars of the Minkowski functionals obtained with 1000 Gaussian simulations. Taken from [461 ]. 



The number of peaks above a threshold v, N(u), can be calculated from the previous differential quantity by integrating 
over the whole parameter space for n and e and over the interval (is, oo) for v. 



N(u) 



V/ ( - 

7T 



1/2 



exp(-z, 2 /2) 7V-1) 



1 — — erfc(7^s) 



+i/ 7 (l - 7 2 )_ exp(- 7 2 i/ 2 s 2 ) + 



(10) 



+t 



1 — -crfc(7^si) 



exp(-7Vi 2 ) 
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where 



= [2(l- 7 2 )p 1/2 , i=(3 - 2 7 2 )" 



1/2 



(11) 



and Nt — (S 1 / 2 *? 2 ) -1 is the total number of peaks on the whole celestial sphere. 

Another interesting quantity is the distribution of the curvature of the peaks. The pdf of peaks with inverse of 
the Gaussian curvature L = k^ 1 between (L, L + dL) above the threshold v, p(L, can be also obtained from 
N(k, e, v) dn de dv: 



P(L) 



1/2 



(2 7 2 ) 4 iL- 5 exp[(2 7 c 2 ) 2 L 



di/exp(-3*V/2)erfc -(2 7 6> 2 L _1 - ^vt 2 



(12) 



The distribution of excentricities can be calculated in a similar manner. The pdf of peaks with eccentricity between 
(e, e + de) above a threshold is, p(e), can be obtained after a straightforward calculation: 



32V6 



1 



(2 



di/exp(z/ 2 /2) |(i/7r) 1 / 2 exp(-G) 1 - edc{H x ' 2 ^vs)/2 

[3ff 2 (l - 7 2 ) 2 + 6i7 3 7 2 (l " 7> 2 + (Hw) 4 ] + 
exp(- S 2 7 V> [5tf 3 7 (l - 7 2 ) 2 ^ + H 4 M 3 (1 - 7 2 )] } , 



(13) 



where H = (2 - e 2 ) 2 /[(3 - 2 7 2 )e 4 + 4(1 - e 2 )] and G = H{ 1 ve 2 ) 2 /{2 - e 2 ) 2 . 

Scalar quantities can be constructed from the derivatives of the CMB field on the sphere. A single scalar can be 
constructed in terms of the ordinary derivative of the field A^. Only two independent scalars can be obtained from 
the second covariant derivatives on the sphere, A ; y. Following [1221 ] . many scalar quantities can be defined from 
the first and second derivatives of the field, associated to the Hessian matrix, the distorsion, the gradient and the 
curvature. However, all except three are correlated. For testing Gaussianity, it is convenient to use normalized scalars 
for which the dependence of the scalars on the power spectrum has been removed. Here, as example, we focus on 
three independent normalized scalars, namely the Laplacian, the fractional anisotropy and the square of th e mo dulus 
of the gradient. The first two scalars have been proved to be very efficient as detectors of non-Gaussianity |l23| . The 
third scalar, the square of the modulus of the gradient, is the only scalar from the list given in that paper which is 
independent from all the others. It depends only on the first derivatives of the field. The Laplacian and the fractional 
anisotropy depend only on second derivatives and can be defined in terms of the eigenvalues Ai, A2, of the negative 
Hessian matrix A of the field: A = (— T ; ij). The eigenvalues, i.e. the negative second derivatives along the two 
principal directions, can be written as a function of the covariant second derivatives of the field A(n): 



Ax = -i 
2 



(A*,W(A*,) -2 (A" 



aw - a ;j ,A ;: 



(14) 



• \ (A : ',) 2(A',A' / A : ',A 



(15) 



We will assume Ai > A2. Considering the values of Ai and A2, three types of points can be distinguished (see e.g. 

[55l ] : hill (both positive), lake (both negative) and saddle (one positive and one negative) 

The normalized Laplacian, or trace of the Hessian matrix, A+, is defined in terms of the eigenvalues as: 



- Ai + A 2 . 

A+ = , —00 < A+ < 00 

er 2 



(16) 



Since the Laplacian is given by linear transformations of the CMB temperature fluctuation field A, if the field is 
Gaussian then its 1-point pdf is also Gaussian: 



P(A+) 



1 

/2~7T 



exp I — — 



(17) 
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The fractional anisotropy [15| has been used in other fields like in the analysis of medical images. The normalized 
quantity f a is defined as: 



fa 



1 A?-A| 



y/2 VA^Ta 

where Ai and A2 are the normalized eigenvalues given by: 



(18) 



Ai 
A 2 



= 1 « ^f^J 



°2 ^/af^l 

The pdf of the normalized fractional anisotropy is given by: 

2/„ 




1 _j_ 1 



P(fa) = 



(l-/a 2 ) 1/2 (l+^) 3/5 



,0</.<l 



(19) 



(20) 



The normalized square of the modulus of the gradient g, which depends only on the first derivatives of the field, is 
defined as: 
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IVAI 



A' l A 



(21) 



where a\ is the dispersion of the unnormalized square of the modulus of the gradient and accounts for the normalization 
factor. In terms of the derivatives of the field with respect to the spherical coordinates (9, <p), g takes the form: 



1 

9 = -5 



dA 



1 



sin 



dA 



(22) 



Taking into account that the square of the gradient modulus is given by the addition of two independent squared 
Gaussian variables, its pdf follows a X2 distribution. Since we consider the normalized quantity g, then its mean and 
dispersion are equal to one and its distribution takes the simple form: 



p(g) = exp(-g) , < D g < 00 



(23) 



The pdfs of the three normalilzed scalars can be seen in Fig. [5] Maps of the same scalars for a random realization 
of a IGRF are shown in Fig. [3l 



III. PHYSICAL EFFECTS PRODUCING DEVIATIONS FROM THE STANDARD IGRF 



There are a number of physical effects which may produce deviations from the standard IGRF. They are normally 
related to secondary anisotropics produced by photon scattering, gravitational effects generated by the non-linear 
evolution of the matter density, variations from standard inflation in the early universe, topological defects, non- 
standard geometry and topology of the universe or primordial magnetic fields. Below we sumarize some relevant 
aspects of the most studied effects. 



A. Secondary anisotropics 

On the last scattering surface CMB anisotropics are generated via the gravitational potential (Sachs- Wolfe effect, 
136]) and the physics of the baryon- photon plasma. These anisotropies are usually refered to as primary anisotropics. 
After matter-radiation decoupling when the temperature drops below w 3000-fC, new anisotropies are generated 
during the trip made by the CMB photons to reach us. These secondary anisotropies can be originated via the 
gravitational redshift suffered by the photons when crossing the gravitational potentials produced by the large-scale 
matter distribution, or by the scattering of the microwave photons with the ionised matter after the reionization 
epoch [z < 10). The gravitational potential can produce two types of effect: the gravitational redshift suff ered by th e 
photons when they cross the varying potential wells formed by the matter evolution (Rees-Sciama effect, |l!5l . Il3l| ). 
and the lensing effect produced by the same gravitational potentials which bends their trayectory (see e.g. Bartelmann 
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and Schneider 2001 for a review). The CMB power spectrum produced by the Rees-Sciama effect is subdominant on 
all a ngular scales [8llll37j . and the expected 3-point correlations have an amplitude much below th e co s mic variance 
[l2ll | . The lensing effect produces a smoothing of the acoustic peaks in the CMB power spectrum [ll7l I X 39t ] . It also 
induces high-order correlations in the CMB temperature and polarization fields [791 ]. 

The scattering of the CMB photons with the ionised medium produces a randomisation of the directions of the 
rescattered photons, implying a supression of the primary CMB anisotropics whose resulting power spectrum depends 
on the Thompson optical depth r as exp(— 2r). In addition, new anisotropics are generated by the scattering with 
the free electrons moving along the line of sight which produces a Doppler effect. This effect is strongly suppressed 



normalised temperature 
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due to the cancellation of velocities along the line of sight. However the Doppler effect can survive cancellation if it 
is modulated by either density or ionization fluctuations. Its amplitude is given by: 



AT 

— = -<tt 



dte- T ^^n e (9,t)v r (6,t) . (24) 



In this formula ot is the Thompson cross-section, n e the electron density and v r the line of sight velocity of the 
electrons. n e can be further expressed as n e (9, t) = n e (t)[l + 5(9, t) + 6i(8, t)] where n e (t) is the mean electron density, 
6 the density fluctuation and Si the ionisation fraction. Two second order effects generating secondary anisotropics 
appear in this formula. The first one is the Dopp ler effect modulated by the density variation and known as Ostriker 
and Vishniac effect in the linear regime [l27l 1 153j ] . The second effect is the Doppler effect modulated by the ionisation 
fraction and is usually referred to as patchy reionization (see e.g. Q). Since they are of second order, both effects 
produce non-Gaussian perturbations in the CMB. However, their amplitude is much smaller than the thermal and 
kinetic Sunyaev-Zeldovich effects (SZ) [l44j generated in the non-linear regime at lower redshifts and produced by 
the scattering of hot, ionised gas associated to collapsed structures. In addition to temperature anisotropies, due to 
the primary quadrupole moment reionisation also produces polarisation. 

For more details about secondary anisotropies the reader is refered to e.g. (sjj and 



B. Non-standard models of the early universe 

1. Non-standard inflationary models 

In the standard, single field, slow roll inflationary model, the dominant linear effects in the evolution of density 
fluctuations preserve the initial Gaussian distribution produced by the quantum fluctuations of the inflaton. On the 
contrary, second-order effects perturb the original Gaussianity of the fluctuations, although resulting in a negligible 
effect in the standard model. In non-standard inflationary models, primordial non-Gaussianity should be added to 
the second-order effects associated to the evolution of density fluctuations after inflation finishes. In any case, it is 
common to characterize phenomenologically the deviations from Gaussianity by introducing the /jvl parameter in 
the gravitational potential [95| : 

<f> = 4>L + fNL<t>L > ( 25 ) 

where <f>i, is the gravitational potential at the linear order and it is distributed following a Gaussian random field. 
In general the non-linear couplig parameter /jvl is a function of the distance vectors and the product is really a 
convolution. However, an effective /jvl which accounts for those complexities can still be used. A detailed study of 
the values of /jvl for several non-standard inflationary models, including multi-field inflation, the curvaton scenario, 
inhomogeneous reheating and the Dirac-Born-Infeld (DBI) inflation inspired by string theory, can be found in [l3| 
(see also [l20j | in this volume). Here it is worth noticing that sometimes the specific /jvl parameter given by eq. 1251 
is called the local non- linear coupling parameter , refering to the fact that is obtained from ^ at the same 
position in space. In addition to /jvl which contains mainly information of the squeezed configurations (those with 
two large and similar wave vectors and the other small), the equilateral non-linear coupling parameter /Jvl* is used 
to caracterize equilateral configurations of the bispectrum for which the lengths of the three wave vectors are equal in 
Fourier space. The two non-linear coupling parameters suppose a fair representation of a large class of models. For 
instance, /j^£ a ' can be generated in curvaton and reheating scenarios whereas /S^" can be produced in DBI inflation 
within the context of String theory. 

In the standard single- field slow-roll inflation the effective /jvl is dominated by second-order gravitational corrections 
leading to values of order unity. These low values require very sensitive measurements to be detected and are not even 
within reach of the Planck mission. Non-standard inflationary models generally predict larger values, some of them 
already constrained by WMAP. A positive detection of /jvl 10 would rule out the standard inflationary models. 

For an ideal experiment with white noise, no foreground residuals and no Galactic and point source mask, it has 
been shown that the optimal estimator for /jvl is based on a bi spectrum test constructed from a cubic combination 
of appropriately filtered temperature and polarization maps |40l [981 Il57l . An extension to this estimator to deal with 
data under realistic experimental conditions has been made by [1581 ] . Several observational constraints on /jvl have 
been derived from different experiments. The first on es w ere obtained with COBE-DMR [3lT[96|. Further constraints 
have been derived with other experiments like VSA [l40j ]. BOOMERANG J^j an d Are heops |46l. l47j . However, the 
best limits have been derived with the various releases of the WMAP data [97l. l99l . Il42| representing an improvement 
of at least an order of magnitud over previous ones (see below) . 
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2. Topological defects 

In standard theories of particle physics the fundamental forces of nature unify progressively when the energy scale 
exceeds certain thresholds. These unified theories imply that the universe went through several phase transitions 
during the early stages of its evolution. At energies above 10 2 GeV, the electromagnetism and the weak nuclear force 
are merged into the electroweak force. At higher energies of the order of 10 15 GeV it is believed that the electroweak 
force unifies with the strong nuclear force, a process usually called Grand Unified Theory (GUT). At even higher 
energies it is speculated that it is possible to merge gravity with the other three interactions. 

Unified theories are based on symmetry. When the symmetry is broken spontaneously via the Higgs mech anism , 
topological defects generically appear [§2] (for a more pedagogical discussion about topological defects see [l49j | ) . 
Depending on the dimensionality of the symmetry which is broken, a type of topological defect is formed. When 
a discrete symmetry is broken domain walls form, in a two-dimensional symmetry breaking cosmic strings appear, 
monopoles form when the symmetry breaking is three-dimensional and textures appear when it is in four or more 
dimensions. In the 1980s topological defects were considered as an alternative scenario to inflationary quantum 
fluctuations for the process of structure formation. CMB observations showed that the former scenario could only 
play a subdominant role as the source of cosmic structure. Moreover, due to the catastrophic effects that the presence 
of domain walls and monopoles would have on our universe only the existence of cosmic strings and textures is usually 
considered. 

A lot of effort has been made to study the cosmological consequences of cosmic strings. It is generally believed that 
the best strategy to detect them is by searching for their imprint on the CMB maps. The non-Gaussian character 
of the density field of strings produce line discontinuities in the CMB anisotropies at arcmin angular scales as a 
consequence of the metric deformation around the strings (the Kaiser-Sttebins effect (89(). However, at larger angular 
scales a Gaussian distribution emerge from the central limit theorem. 

Cosmic textures were first studied in detail by [l47j ]. They are much more diffuse than the other defects which are 
localized at a point (monopolc), on a line (cosmic string) or a surface (domain wall). Contrary to the others, they are 
unstable and consist of twisted configurations of fields which collapse and unwind. Each texture creates a time- varying 
gravitational potential which produces red- or blue-shift to the CMB photons passing through such a region. Thus, 
textures generate hot and cold spots on the CMB anisotropy maps whose amplitude is set by the symmetry breaking 
energy scale. The sh ape o f the spots is approximately spherically symmetric and an approximated analytical formula 
has been derived by |l48| . Recently, a very cold spot detected in the WMAP temperature map has been found to be 
consistent with the effect produced by a texture [44( . If confirmed, this result will have outstanding consequences in 
our understanding of the universe (see Sec. I VI F~2l for more details on this finding). 

C. Non-standard geometry and topology 

The geometry of the observable universe is believed to be well approximated on large scales by a Friedmann- 
Lemaitre-Robertson- Walker (FLRW) metric characterized by a homogeneous and isotropic space-time. Support for 
the homogeneity of the universe is provided by the largest surveys of the galaxy distribution (e.g. SDSS DR6 [H, 
2dFGRS |3y|) and by the smallness of the CMB temperature fluctuations. The cosmological principle, stating that 
the universe is homogeneous and isotropic on large scales, follows from those observations and the assumption that 
we are not located at a special place in the universe (the Copernican principle). 

However, recent observations of the CMB as measured by the WMAP satellite might question the cosmological 
principle. They are the large scale power asymmetry found between the two hemispher es of a reference frame close 
to that of the ecliptic one [53, l72 | . t he pl anarity and alignment of the low multipoles [49|, Il04l | , the non -Gaussian cold 
spot present at scales of sw 10° |42l. I X 50j ] and the alignment of local structures also at similar scales [l5l[ Il55j (see 
Sec. IVII for more details). These results may imply the existence of priviledged directions in the CMB map and thus 
motivates the study of alternatives to the standard FLRW metric. 

An interesting class of alternative models is that for which the metric is homogeneous but anisotropic. They are 
known as Bianchi models and are classified acording to their space-time properties. Their predictions for the CMB 
anisotropy were studied in [Io| . Since the signatures left by those m odels appear on large angular scales, they have been 
already constrained by the COBE-DMR experiment data [13, [sH, Ill6| | . One particularly interesting case is Bianchi 
Vllh which experiences anisotropic expansion and global rotation. These two properties produce a characteristic 
pattern in the CMB in the form of a spiral pattern an d sp ots. This model has been recently used to account for some 
of the large scale WMAP anomalies mentioned above [24], HH, HE H3 (see Sec. lVIFTI for more details). An important 
problem with this model comparison approach is that the CMB anisotropies are not computed in a self-consistent 
way but are assumed to be the sum of two independent components: an isotropic one produced by the energy-density 
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fluctuations which is assumed to be the same than for the FLRW model, and an anisotropic component which is the 
determinisitc effect produced by the anisotropic model. 

A different source of anisotropic features in the CMB sky is the global topology of the universe. The local character 
of the General Theory of Relativity does not theoretically constrain it. If the topology of the universe is non-trivial 
(i.e. multiconnected, meaning that there is not a unique way to connect two points by geodesies) then CMB photons 
originated from the same location on the last scattering surface can be observed in different directions. This effect 
manifests itself in the CMB sky as anisotropic patterns, correlated (matched) circles [3!| or, more generically, as 
devi ation s from a IGRF (for a discussion on different topologies and tests developed to detect them see the reviews by 
[l03j and [lOS]) . An additional consequence of a multiconnected universe is the lack of fluctuations in the CMB above 
the wavelength corresponding to the size of the universe. This property has led to several authors to suggest that the 
low quad rupole a nd the alignment of the low multipoles measured by WMAP might be an evidence o f a non-tr ivial 
topology [M[m|]. Constraints on the to polo gy of the universe started with the COBE-DMR data [H, [HI [l33j] and 
followed with the WMAP data [H HoU Il02l OH, EU- All those analyses concluded that the WMAP data do not 
show any evidence of multiconnected universes. 

D. Primordial magnetic fields 

In recent years there has been an increasing interest in studying the consequences that the possible existence of 
primordial magnetic fields might have on the CMB. Magnetic fields of order of a few jiG have been measured in a 
wide range of astrophysical structures, from individual galaxies [71| to galaxy clusters [lOOj ]. It is also widely beleived 
that they are present in superclusters. Although the origin of those magnetic fields is still unclear and their existence 
does not necessarily imply a primordial origin, studying the interplay between magnetic fields and CMB is justified 
by the important consequences that it may have on the CMB temperature and polarization anisotropy, and also on 
the distorsion of the blackbody spectrum (for reviews on this topic the reader is referred to (58l. loTj). In particular, 
the presence of magnetic fields introduce non-Gaussianity in the CMB anisotropy since its amplitude depends on the 
square of the magnetic field intensity. 

Two very different cases can be considered for the primordial magnetic field: a uniform field and a stochastic 
one. The former breaks the spatial isotropy of the background geometry by introducin g sh ear through an anisotropic 
stress. This leads to the well known homogeneous and anisotropic Bianchi models [lOj. Since those models also 
imply an anisotropic CMB field, a uniform magnetic field generates phase correlations between different ai m . The 
uniform magnetic field has been constrained with the CMB quadrupole and also with the phase correlations, implying 
comparable constraints on the magnetic field intensity of a few nG [llj, |33|, |57| . 

The stochastic magnetic field is a more realistic scenario which can be generated during inflation. In this case the 
isotropy of the background geometry is preserved. Temperature and polarization anisotropics are generated through 
the vector and tensor modes associated to the magnetic field energy-momentum tensor. Allowed amplitudes of the 
magnetic field intensity of about several nG can produce a potentially observable B-mode polarization for nearly scale 
invariant spectra [1091 ]. This signal could be distinguished from the one generated by the inflationary gravitational 
wave background by its non-Gaussian character. In addition Faraday rotation induces B-mode polarization from the 
ordinary E-mode with the characteristic v~~ 2 dependence. 

Therefore, the presence of a primordial magnetic field can leave unambiguous inprints in the CMB anisotropy that 
would allow its identification with the sensitive data expected from the coming experiments. 

IV. METHODS TO TEST GAUSSIANITY 

Testing the Gaussianity of CMB data is not an easy task. In principle it consists in just proving the properties of 
the IGRF that we discussed in Sec. Ill Bl isotropy and multinormality. The CMB data represent a single realization of 
the underlying random field which for the standard model is nearly Gaussian and isotropic. In practice the analysis is 
complicated by the characteristics of the experiment which need to be known very precisely: calibration uncertainties, 
instrumental noise (white and 1/f) which is normally anisotropic in pixel space depending on the scanning strategy, 
beam response (usually close to Gaussian), data processing, ... And by foreground contamination which demands 
certain previous cleaning operations in the data, requiring masking certain areas where the foregrounds are very 
intense and leaving some amount of residuals in the rest of the surveyed area. The result of cleaning depends on our 
a priori knowledge of the physical properties of the foregrounds and the component separation method used for their 
removal. These ingredients must be considered in the analysis by performing simulations accounting for them. 

Unless one is interested in the compatibility of the data with a specific alternative model for which and optimal 
method may be found (as specific non-standard inflation, geometry or topology), there are infinite ways in which a 
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random field can deviate from the IGRF one. Methods to test Gaussianity can be classified by the type of property 
that they try to prove. Typical examples are cumulants and n-point correlation functions in real space (the former 
should vanish and the latter either vanish for the odd order or can be expressed in terms of 2-point correlation functions 
for the even order), moments in spherical harmonics space (bispectrum, trispectrum, ...), or moments in other spaces 
to which the data is transformed by linear operations which preserve Gaussianity: filters, wavelets, signal-to-noise 
eigenvectors, ... Other approaches may test different properties of the CMB random field, like the morphology of 
the data using the Minkowski Functionals or the geometry of the data using scalar quantities constructed by the 
first and second covariant deriv ative s as, for example, the local curvature (see Sec. Ill Bp . These and other methods 
are described in more detail in [l!9l | and Q and references therein. As example of the typical statistical procedure 
followed in the analysis, below we focus on a few statistical methods which have been often used in the literature 
based on the Minkowski Functionals, bispectrum or wavelets. 

As we have described in Sec. Ill Bt for the excursion set above a given threshold v, there are t hree Minkowski 
Functionals on the sphere: total contour lenght C(v), total area A(v) and the genus G(v) (see [l38| ). In the case of 
an IGRF their expected values follow simple analytical expresions as a function of v (eqs. l4l5t6|) . However, simulations 
are needed to account for the experimental characteristics, basically noise, beam response and mask. As it can be 
shown with simulations, the 1-pdf for any Minkowski Functional at each v follows a nice bell-shape distribution, 
implying the natural choice of a generalized x 2 as the appropriate statistical test to be used in this case to combine 
all the information. More specifically, considering n t h different thresholds we can define a 3n t h vector v, 

v = {A{v x ), .., A{v nth ), Cfa), .., C(u nth ), G{v x ), .., G(v nth )) . (26) 

The generalized x 2 statistic to test the Gaussianity of a data map can be then constructed as 

3nth 
i J = l 

where <> is the expected value for the Gaussian case and CV, is the covariance matrix, Cy =< ViVj > — < V{ >< v j >, 
both of them usually constructed with simulations. For testing the compatibility of the data with a non-Gaussian 
model (e.g. a non-standard inflationary model characterized by the Jnl parameter) we simply have to use the 
corresponding expected value and covariance matrix for the Minkowski Functionals at different thresholds in eq. 1271 
If the deviations from Gaussianity are small (as in the case of the Jnl models with /a^ < 1000) the covariance matrix 
can be well aproximated by the one of the Gaussian case. More information is added to the analysis by considering 
yires different resolutions of a given data map. Including this extra information simply increases the vectors and 
covariance matrix present in the x 2 expression to a dimension 3n t ^n res . Examples of applications of this method to 
different data sets can be found in \91§ for the 1-year WMAP data, [H[ for the BOOMERanG 3-year data and H3 
for the Archeops data. Recently, perturvative formulae of the Minkowski Functionals as a function of Jnl have been 
derived for the /jv l models 174] . The results of applying them to the 3-year WMAP data show constraints on Jml 
very similar to the ones obtained with the optimal bispectrum 75] . 

Generically, non-standard models of inflation produce small deviations of Gaussianity which are more promi- 
nent in the 3-point correlation function or equivalently, its harmonic transform the bispectrum B™f™* m3 =< 
o-i imi O'i2m2 a hm3 >• The trispectrum, the harmonic transform of the 4-point correlation function, can also play 
an important role in discriminating inflation models since some models do not produce any bispectra but produce 
significant trispectra, or produce similar amplitudes of the bispectra but very different trispectra (e.g. DBI inflation, 
[78j |. or New Ekpyrotic Cosmology, [I(|). Here we briefly describe the bispectrum, for more details on it and the 
trispectrum see |94j |. The average bispectrum, i?; 1 / 2 ; 3 , is the rotationally invariant third order moment of spherical 
harmonic coefficients and is given by the following expresion (see [80j and [l3j for more details on this and higher 
order moments): 

B hhh = X! (mi m 2 m 3 ) < a hm 1 ai 2m2 ^hm 3 > , (28) 

mim2m3 ^ ' 

where (...) is the Wigner-3j symbol. The bispectrum must satisfy the selection rules. Rotational invariance of the 
3-point correlation function implies that the bispectrum can be written as 

Kzr ni = G z\:ir^H , (29) 

where G™^™ 2 ™ 3 is the Gaunt factor and &; 1 i 2 i 3 is a real symmetric function of l\, l 2 and I3 called the reduced 
bispectrum. By substituting eq. [2H in eq. [55] it is straight forward to obtain the following relation between the 
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averaged bispectrum and the reduced bispectrum 



(2/ 1 + l)(2/ 2 + l)(2/ 3 + l) (l x l 2 l 3 



4tt 



Khh ■ (30) 



Therefore under rotational invariance the reduced bispectrum contains all physical information of the bispectrum. In 
particular &i 1 / 2 ; 3 is fully identified by the inflationary models characterized by the non-linear coupling parameter /jvx- 
Thus these models can be optimally tested by using the averaged bispectrum. The computation of the full averaged 
bispectrum scales as N^, where N pix is the number of pixels, and is already not feasible for the WMAP data. [98| 
solved this problem by constructing a cubic statistic, the KSW estimator, that combines the triangle configurations of 
the bispectrum optimally for determining Jj^ff' , with its computation scaling as . The extension of this estimator 



to /jvi"^ and also including polarization, has been made in [40j | and 157 1. The strongest constraints on the non- linear 



coupling parameter up to date have been obtained by applying the KSW estimator to the WMAP data (see next 
section). 

Wavelets are compensated filters which allow one to extract information which is localized in both real and harmonics 
space. In particular wavelets may provide information on the position and scale of different features in astrophysical 
images. They can be more sensitive than classical methods (for a review on wavelets see [88[ and for applications to 
the CMB see [152]). One example which has been used many times in cosmological applications is the Mexican hat 
wavelet defined as the Laplacian of a Gaussian function. For CMB analyses the extension of an Euclidean wavelet 
to the sphere can be made with an inverse stereographic projecti on [g| . By this procedure the spherical Mexican hat 
wavelet can be constructed preserving dilation and compensation [l 18j ] . and was first applied to CMB analyses by (30j 
using COBE-DMR data. Also it can be made directional, i.e. sensitive not only to the scale but also to the orientation 
of a feature, by simply considering different widths along the two axes of the original 2-dimensional Gaussian. 
A very useful property to study directional properties of structures in an image is the so called steerability. In general, 
a directional or non-axisymmetric filter is said to be steerable if any rotation about itself can be expressed as a finite 



linear combination of non-rotated basis filters. This concept has been recently extended to the sphere by 154]. An 
important consequence of the steerability property is that it makes the wavelet analysis very efficient computationally 
(see Sec. lVIDj) . An interesting spherical steerable wavelet is the second Gaussian derivative one which may be rotated 
in terms of three basis wavelets: the second derivative in direction x, the second derivative in direction y, and the 
cross-derivative. This wavelet provides information on the three local morphological measures of orientation, signed 
intensity (amplitude at the orientation which maximizes the absolute value of the coefficient) and elongation. It has 
been recently applied to the CMB analysis to test global isotropy and Gaussianity jl5ll . Il55l Il56l ] (see Sec. lVIDj) . 



V. CONSTRAINTS FROM OBSERVATIONS 



Previously to any Gaussianity study, one major problem in the analysis of the observational data is the separation 
of the different Galactic and extragalactic components from the CMB itself. A key property to distinguish them 
is the frequency dependence of the specific emission in the microwave range. Thus, the intensity of the Galactic 
synchrotron and free-free emissions decreases with frequency as a power-law with approximated spectral indexes 3 
and 2, respectively. On the contrary, the intensity of the Galactic thermal dust increases with frequency following 
approximately a grey-body spectrum with emissivity oc v 2 and temperature Tu ~ 10 — 20K. Extragalactic sources 
emitting in the microwave band are typically radio galaxies and IR galaxies, with their emission dominating at 
frequencies < 100GHz for the former (synchrotron-like) and at higher ones (dust-like) for the latter. An additional 
extragalactic emission comes from galaxy clusters through the SZ effect, as was discussed in Sec. MI Al The spatial 
distribution of the different foregrounds is also very different from the CMB one, differing very much from that of 
an IGRF. Thus, the Galactic foreground emissions dominate at large angular scales with the power spectrum of 
fluctuations decaying approximately as Ce oc £~ 3 . On the contrary the extragalactic foregrounds dominate at the 
smallest angular scales and appear as point sour ces for typical CMB experiments. (For more details on the properties 
of foregrounds see contributions by [53| and [l29t ] to this volume) . All these properties have to be exploited in order to 
best disentangle the foreground emissions from the CMB one. Indeed, the numerous component separation methods 
already developed take advantage of our knowledge of those foreground properties to obtain a clean CMB map (see 
[511 and [j| for further reading on this topic). 

The Gaussianity of the CMB signal has been studied with data measured by many different experiments. The 
first systematic analysis was carried out with the all-sky COBE-DMR satellite [la , 1 14lj | . Other analyses were based 
on data covering a fraction of the sky and obtained from experiments onboard stratospheric balloons, like QMAP, 
MAXIMA, BOOMERanG or Archeops, and from ground-based ones like Saskatoon, QMASK or VSA. The result 
of most of those Gaussianity studies was a systematic compatibility with the standard IGRF. Deviations were also 
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claimed in a few cases which were later proved to be due to either systematics or an incomplete analysis. Upper limits 
on the fx L al para mete r were derived of approximately a few thousands [U, [9(| . For a more detailed discussion on 
those analyses see jl!9j . A significant improvement has been recently achieved with measurements by BOOMERanG 
and Archeops, lowering the upper limit to Jnl ^5 1000 at the 2a c.l. [I?], [H|. It is worth mentioning a deviation 
from the IGRF found in one of the VSA fields, the Corona Borealis supercluster [65|, Il34j . The deviation consists 
in a strong and resolved negative spot, of ~ — 250/iK and angular size of ~ 20 arcmin, which is not associated with 
any of the clusters of that supercluster. A SZ effect produced by a diffuse, extended warm/hot gas distribution has 
been suggested as a possible explanation [66]. This hypothesis, if confirmed, would be of relevance for providing the 
location of the missing baryons in the local universe. 

The precision with which the Gaussianity of the CMB can be tested has been strongly improved with the quality data 
measured with the WMAP satellite [U [76|, [z3| ■ The WMAP team has provided "cleaned" CMB maps at the Q, V and 
W frequency bands for the 5 year data collected by the experiment, and masks covering a region around the Galactic 
plane and a catalogue of radio sources where the foreground emission cannot be removed at the required sensitivity. 
By a noise- weighting combination of the cleaned maps at the cosmological frequencies and applying a conservative 
mask the WMAP team has performed a Gaussianity study of the data based on the Minkowski Functionals and the 
optimal bispectrum. F or bo th quantities the 1, 3 and 5-year WMAP data releases have been found to be compatible 
with the IGRF [9?], IH, Il42( . Stringent limits have been also derived on the local and equilateral f^L parameter using 
an optimal bispectrum- like quantity: -9 < f§ c L al < 111 and -151 < Ql < 253 at the 2a c.l. 



VI. WMAP ANOMALIES 



The WMAP team found the WMAP data consistent with the IGRF using the Minkowski functionals and the 
bispectrum. Subsequently, many works have tested the Gaussianity of the WMAP data in many different ways. Some 
of them have found agreement with Gaussianity whereas others found significant deviations of the IGRF. Examples of 
the former are analyses based on 3-point correlation analysis [341. |64| . i ntegrated bispectrum [28[ , real space statistics 
[Tij or isotropy analyses based on the Bipolar power spectru m |69L|70H. Examp les o f the latter are analyses based on 
phase correlations [3al . the genus |l28j , isotropic wavelets [HI, Il25l . Il50t| , 1-pdf |l24j , isotropy an alyses based on loca l 
n-point correlations [59l| , local curvature [73[ , multipole vectors |37l |38| or directional wavelets [l!2l 1 15 lL Il55l | l56l | . 
However, some of the analyses have been performed on the whole sky internal linear combination (ILC) map [1 81 ] or 
similar maps (e.g. [l45j ]) which is well known to suffer from Galactic contamination and, as stated by the WMAP 
team itself, should not be used for cosmological analyses. From here on, we will concentrate on the most relevant 
works based on WMAP maps where a certain region around the Galactic plane has been masked as well as several 
hundred of extragalactic sources. The deviations or "anomalies" reported have been detected using other statistical 
quantities different from the ones or igina lly used by the WMAP team . They are the North-South asymmetry p59l.l72l]. 
alignment of the lo w multipoles [49l Il04| . the cold spot [42l. Eil [l5"oT | . non- G auss ian features detected with directional 
wavelets [112| . 1 1 13| , alignment of CMB structures 15lL 155j . low variance jl24| . Below we describe these anomalies 



and discuss several relevant aspects like significance, origin, etc. 



A. North-south asymmetry 



North-South asymmetries in ecliptic coordinates have been observed in the WMAP CMB maps using several 
local quantities: power spectrum and 2 and 3 -point correlation functions [59L [sTJ , Minkowski functionals [60j , local 
power spectrum [54]> [zl], local bispectra [l05j ] and local curvature [73j |. Varying the coordinate system, a maximum 
asymmetry for the two hemispheres is obtained for a system whose north pole is lying at (9, <f>) — (80°, 57°). This 
direction is close to the north ecliptic pole (9, </>) = (60°, 96°) (see Fig. 2]). The result of all those works basically 
indicates a significant lack of power in the north ecliptic hemisphere compared to the south one. The asymmetry has 
been also confirmed by fl9l . [20 ] using a pair angular separation histogram method. 

The asymmetry remains stable with respect to variations in the Galaxy cut and to the frequency band. Also a 
similar asymmetry is found in the COBE-DMR map with the axis of maximum asymmetry close to the one found in 
the WMAP data. Analyses of the possible foreground contamination and known systematics do not seem to be the 
cause of such asymmetry (59[. 

More recently, the asymmetry has been found again in the 3-year WMAP CMB map [62j . 




FIG. 4: Directions in the microwave sky derived from the anomalies found in WMAP data: northern direction of the North-South asymmetry (NSA), axis of evil 
(AoE), the cold spot (CS), cluster of positive total weights (CpTW), perpendicular axis to the positive total weights plane (pTW), perpendicular axis to the negative 
total weights plane (nTW). For reference the North ecliptic pole (NEP) and the northern direction of the CMB dipole (NDD) are also shown. 
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B. Alignment of the low multipoles 



The lowest mul poles. especially I — 2, 3, of the WMAP data have been found to be anomalously planar and aligned 
[22l [37l l38l l49l. [90l Il43j . In Fig. [5] the quadrupole and octopole of the ILC map are shown. Both multipoles present 
maxima and minima following a planar shape, whose perpendicular axis points towards similar direction called axis 
of evil. The axes of the two multipoles are separated by m 10°. The probability that the two directions are separated 
by t hat angle or less by chance is « 1.5%. Further alignments have also been claimed for higher multipoles I < 5 
[l04| and I = 6, 7 (63|. The northern end of the alignment points towards (9, <fi) — (30°, 260°), a direction close to the 
CMB dipole one whose northern end is at (0, </>) = (42°, 264°) (see Fig. |J). 

A problem which appears when trying to estimate the low multipole components is that they are very much 
affected by the mask. Varying the mask produces significant changes in their amplitude estimates, especially for the 
quadrupole, implying consequently uncertainties in the determination of their axes. Detailed analyses of this effect 
tend to weaken the significance of the detection Il07l | . 



C. The cold spot 



A large and prominent cold spot has been observed in the WMAP data which is hard to explain within the standard 
inflationary scenario. It was detected with the SMHW as defined in [l!8j |. The first evidence came from the kurtosis 
of the wavelet coefficients of the first-year data which showed an ex cess w ith respect to the Gaussian hypothesis at 
a wavelet scale of ~ 4° (corresponding to a structure size of « io°) a very cold spot at Galactic coordinates 

(9, </>) = (—57°, 209°) was identified as the possible source of the excess (see Fig. [6]). An analysis of the area of the 
spots at different thresholds in the SMHW coefficient map at around 4° proved that indeed the cold spot had a very 
large area and was the source of the excess of the kurtosis [42[ . 

The cold spot has been also shown to deviate from Gaussianity using the Max and Higher Criticism estimators 
[32I l42l ]. A s tudy of the morphology of the spot with the elliptical MHW on the sphere has found an almost 
circular shape [43] ■ In the same paper the possible foreground contribution was considered in detail, concluding that 
contributions from the SZ effect and the Galaxy had to be negligible. 

The cold spot has been also identified as the most prominent spot in the CMB sky using steerable wavelets [151j . In 
that work two other spots are identified as deviations from the IGRF three- year WMAP best fitting model. However, 
the deviation from Gaussianity seen in the kurtosis of the coefficients of that wavelet [l56j | cannot be assigned 
ex clusively t o those three spots. 

[112L Ill3l ] have also detected a number of non-Gaussian features, including the cold spot, using the directional 
wavelets elliptical Mexican hat and Morlet. 

All the previous results have been confirmed with the 3-year WMAP CMB map [3| and are expected to be almost 
unaltered for the 5-year data. 



D. Alignment and signed-intensity of local structures 

Some of the previous anomalies imply preferred directions in the sky that, under the assumption of Gaussianity, 
represent deviations from statistical isotropy. Here we describe a different violation of statistical isotropy based on 
the alignment of CMB structures. The structure s are identified by convolving the CMB map with the steerable 
wavelet formed by the second Gaussian derivative [154| . For each scale and position in the sky the wavelet identifies 
the orientation which maximizes the absolute value of the wavelet coefficients. Thus this orientation corresponds to 
the characteristic orientation of the l ocal feature of the signal. The wavelet coefficient in that specific orientation 



defines the so-called signed-intensity [15l| . It should be remarked that these two quantities, local orientation and 
signed-intensity, are computationally feasible because of the steerability property (see Sec. IIV[) . 

Once the local orientation is determined for every position in the sky at a given scale, we can construct the following 
isotropy test. First, the great circle passing by that position and tangent to the local orientation is identified. Every 
other pixel which is now crossed by that great circle is considered to be seen by the local orientation of the first position 
(this procedure is illustrated in Fig. [7]) , with a weight naturally given by the absolute value of the signed-intensity of 
the first position. Second, we can assign to each pixel in the sky the total weight given by the sum of the weights 
of all the pixels which see that pixel. This new total weight signal is even under parity and thus its analysis can be 
restricted to one hemisphere of reference. The highest total weights represent the positions towards which the CMB 
features are predominantly directed while the lowest ones represent the positions predominantly avoided by the CMB 
features. Of course, for an all-sky IGRF all the pixels should have the same total weight on average. 
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Th is isotropy test based on the alignment of local orientations has been applied to the first [155j | and three-year 
[l5ll ] WMAP data. In both cases a significant violation of isotropy was found at a scale around 10°. The highest 
total weights (above 3c) define an axis located very close to the ecliptic one. The highest and the lowest total weights 
define two planes whose normal axes are close to the CMB dipole one (Fig [5]) . 

Besides the alignment of the local orientations, an analysis of the signed-intensities show three spots at angular 
scales also around 10° containing a total of 39 1°. 8-pixels whose values have a (signal+noise) probability in t he 3tr 
tails (see Fig. This pattern is similar to the one found with the axisymmetric Mexican hat wavelet jl5dl |. The 
three spots are located in the southern galactic hemisphere confirming the north-south asymmetry. Two of the them 
are cold, one being identified with the cold spot already detected in jl5dl |. 



The two anomalies which appear at similar scales, with global significance levels around 1 per cent, are however 
quite independent. 



Very recently [124j | the 1-pdf of the WMAP three-year data has been analysed finding an anomalous low value 
of the variance as compared to the one expected from the WMAP best-fit cosmological model. The result is even 
more prominent if only the north ecliptic hemisphere is considered (see Fig. I10[) . in agreement with the lack of power 
found in that hemisphere by previous works (see e.g. The variance of the CMB signal is obtained by fitting the 

normalized temperature distribution to a Gaussian of zero mean and unit variance. The significance of the result is 
around 1% (see Fig. [TTj) . 

In order to find a possible origin for this anomaly the behaviour of single radiometer and single year data as well 
as the effect of residual foregrounds and 1/f noise have been studied. None of these possibilities can explain the low 
value of the variance. 

Since the largest contribution to the variance comes from the lower multipoles, it is interesting to see if the low 
quadrupole measured by COBE and WMAP is the cause of the anomalously low variance. Performing the same 
analysis after subtracting the best-fit quadrupole outside the KpO mask the significant of the result is slightly reduced 
although the variance is still anomalously low. 

Beyond the inconsistency found between the best-fit model and the measured variance, one could ask if the latter 
is consistent with the actual measured power spectrum of the WMAP data. The analysis performed by the same 
authors show that a strong discrepancy is indeed found. This last result suggests a possible deviation of the CMB 
data form the IGRF. 



Given that neither foreground contamination nor known systematics seem to be causing most of the previous 
anomalies, there has been a number of attempts to explain the cause of the W MAP ano malies by an intrinsic origin. 
Among them we mention the Rees-Sciama effect produced by large voids [H, Il35l Il46l | , inhomogeneous @, Il06l | or 
anisotropic universes [85| and cosmic defects [44j . Although no further evidence has been found for those explanations 
most of them still remain as plausible. Below we discuss in more detail two interesting possibilities: the Bianchi model 
and the cosmic texture. 



A first interesting attempt to explain the best studied WMAP anomalies was performed by [85j, |86j who suggested 
that such features might be produced by an anisotropic universe, the Bianchi VII^ model. This type of model 
has a global anisotropic expansion and vorticity that produces geodesic focusing and a spiral pattern in the CMB 
anisotropy at large angular scales. By fitting the free parameters of that model to the WMAP map the large scale 
CMB anisotropies produced by its non-standard geometry are determined. After subtracting that pattern from the 
WMAP data, the best studied WMAP anomalies, namely the low-multipolc alignments, the north-south asymmetry 
and the cold spot, were significantly reduced. 

This result seemed to suggest that the universe was not well represented by the homogeneous and isotropic FLRW 
model and that a perturbation in the form of a Bianchi VII?! was a better representation. However, a more detailed 
examination of the best-fitted parameters of the Bianchi model when a dark energy was included showed values of 
the dark energy and matter energy density far from the ones measured by many current cosmological tests [24l |87| . 




E. Low variance 



F. Cosmological consequences 



1. Bianchi model 
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Considering other Bianchi models with vorticity and shear, like the Bianchi IX with a closed geometry, do not help 
since they do not exhibit geodesic focusing or the spiral pattern. 

Very recently, f2f| computed the Bayesian evidence of the Bianchi template when the cold spot was not included in 
the analysis (see below for an alternative interpretation of the cold spot). The result was that the evidence was now 
significantly reduced, reinforcing the idea that the cold spot was likely to be driving any Bianchi VII;, detection. 



2. The cold spot texture 

Recently, [3] performed a Bayesian evidence analysis to test the hypothesis that the cold spot is produced by 
a cosmic texture. This is a type of toplogical defect that, as explained in Sec. IIII B 2| in the standard theories of 
unification of the fundamental forces of nature, are expected to appear in the early universe. The motivation for 
considering a texture as the origin of the cold spot is its typical spherical anisotropy pattern left in the CMB and the 
relatively small number of spots expected at scales of several degrees (see Fig. [T"2|) . 

The hypothesis test considered consisted in the comparison of the following two hypotheses: the null hypothesis Hq 
for which the cold spot is just a rare fluctuation of the IGRF predicted by the standard inflationary model, and the 
alternative hypothesis, Hi for which on top of the inflationary Gaussian CMB fluctuations there are non-Guassian 
ones as produced by the texture model. The test is applied to a circular area of 40° diameter centered on the cold 
spot position. The posterior probability ratio of the two hypotheses is: 

PKgijD) Eh Prjm) 
9 Pr(H \B) E Pr(Ho) ' [ ' 

where D is the data vector, and Ei the evidence which is the average likelihood L with respect to the priors 11(0^) in 
the parameters 0^ of the hypothesis Hi, 

Ei = Pr{Ti\Hi) = J L i (6 i |ff i )n(e i )d9i , (32) 

and Pr(Hi) is the a priori probability of hypothesis Hi. The a priori probability ratio is usually set to unit for lack 
of information, but since in our case the analysis is centered in the cold spot (an a posteriori selected pixel) it should 
be given by the sky fraction f s covered by textures. Given that a scale-invariant distribution of spots is expected and 
considering only textures above 1° (photon difusion would smear out textures smaller than that) f s — 0.017. The 
likelihood funtion is simply L cx (~x 2 / 2 ) where \ 2 = ( D - T)'N _1 (D - T) and N is the CMB+noise covariance 
matrix. The anisotropy pattern T produced by a texture can be approximated by an analytical spherical profile |l49l | 
with only two free parameters: the amplitude which is related to the fundamental symmetry-breaking energy scale, 
and the angular scale which is related to the redshift when the texture unwinds. 

The result of the analysis is a probability ratio p — 2.5 for the three-year WMAP data, favouring the texture 
hypothesis (see Fig [T3] for the resulting best-fit texture template). This result is slightly increased (p = 2.7) for the 
five-year data due to the reduction in the noise amplitude. The texture interpretation helps to alleviate the excess 
found in the kurtosis, being at the same time compatible with the observed abundance, shape, size and amplitude of 
the spot. In particular, the symmetry-breaking scale inferred from this analysis, 4>q m 8.7 x 10 15 GeV, is compatible 
with the upper limit obtained from the CMB power spectrum analysis |2lj . 

More recently, [!|| have extended the Bayesian evidence analysis to models based on the Rees-Sciama effect produced 
by voids or on the SZ effect. The result is that, contrary to the texture model, no positive evidence is found for those 
models. 



VII. CONCLUDING REMARKS AND FUTURE PERSPECTIVES 



The main conclusion derived from the large amount of analyses, involving a variety of methods, performed to test 
the Gaussianity of the CMB anisotropies, specially those measured by WMAP, is that the standard IGRF prediction 
is a good representation of their properties as a first approach. Furthermore, a number of significant deviations from 
the ideal IGRF has also been reported whose origin and interpretation is still under debate. Some of them might have 
to do with foreground residuals or unknown systematics while others could be a hint of new physics with profound 
implications for our understanding of the universe. Probably, to answer those questions we will have to wait for new 
data coming from the advanced experiments being built and expected to be operative in the next years. 

Future experiments will shed light on the open questions remaining from the up-to-date analyses of the CMB data, 
specially on the WMAP anomalies discussed in Sec. I VII In particular, the Planck mission is expected to provide 
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all-sky, high quality, multifrequency maps in the frequency range 30 — 900GHz. The wider frequency range and 
the higher sensitivity and resolution will allow an improvement in the quality of the resulting CMB map. As a 
consequence, an improvement in the control of the foreground emission is expected as well as a reduction in the sky 
area required to be masked. In addition to the temperature, improvements are also expected in the polarization maps 
which will be provided by Planck, meaning an important complement for probing the nature of the anomalies as well 
as for testing the different physical interpretations proposed for them. Missions for measuring polarization at the 
highest sensitivity allowed by present technology, and with the main aim of probing the existence of the gravitational 
wave background, have been recently proposed to both agencies, ESA and NASA. As for the temperature, the linear 
polarization expected to be produced as a consequence of the standard inflationary period of the universe also posseses 
properties very close to those of the IGRF studied in Sec. [TTJ Therefore, extensions of the methods already discussed 
in Sec IIVI for temperature are also expected to be applied to test the Gausianity of the future polarization maps. 
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FIG. 6: The cold spot in the 5-years WMAP combined map. This map is a noise- weighted combination of the V and W maps 
given in [77| where the pixels contaminated by Galactic or extragalactic foregrounds have been masked with the KQ75 mask. 
The WMAP maps and masks can be found in the LAMBDA web page. 
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FIG. 7: A simulated map of the signed-intensity where the great circles corresponding to two local features are shown. All the 
pixels crossed by each great circle are considered to be seen by the corresponding local feature. The position which is crossed 
by both great circles is thus said to be seen twice. See [l5l| for more details. 
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FIG. 8: The highest (lowest) total weights above (below) 3cr are plotted in red (blue). A cluster of positive total weights can 
be seen close to the NEP. The red (blue) dashed line is the best fit plane to the highest (lowest) total weigh ts. T he normal 
axes defined by both planes point towards directions close to the northern direction of the dipole NDD. From [l5ll ]. 




FIG. 9: Map of signed-intensities with probabiliy in the 3<r tails. In blue are the pixels with negative values and in red the 
positive ones. The signed-intensities are grouped around three clusters, one formed by the negative values and other two by 
the positive o nes. Also shown are the positions of the South ecliptic pole (SEP) and the southern direction of the CMB dipole 
(SDD). From [l5ll ]. 
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FIG. 11: Theoretical pdf of the CMB variance (calculated following [2j|, solid line) compared to the averaged pdf obtained 
from 60000 Gaussian simulations of the WMAP best-fit model over the whole sky (blue dot-dashed line) and using only the 
pixels out side the KpO mask (red dashed line). The solid vertical line indicates the value obtained from the 3- year WMAP 
data. See 124] for more details. 



28 




FIG. 12: Image of the high resolution texture simulation performed by N. Turok and V. Travieso and available at 
(http://www.damtp.cam.ac.uk/cosmos/viz/movies/neil.html). The distribution of texture spots on the sky is predicted to 
be scale-invariant. 
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FIG. 13: A 43° x 43° patch centered at Galactic coordinates (—57°, 209°) and obtained from the 5-year WMAP combined map 
given in Fig. [6] is shown in the top panel. The best fit texture template is in the middle and the result of subtracting it from 
the WMAP map in the bottom. The units shown in the colorbars are rriK. See [3] for more details. 



